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Abstract 

This work is devoted to the experimental study of the longitu- 
dinal hadronic shower development in the ATLAS barrel combined 
prototype calorimeter consisting of the lead-liquid argon electro- 
magnetic part and the iron-scintillator hadronic part. The results 
have been obtained on the basis of the 1996 combined test beam 
data which have been taken on the H8 beam of the CERN SPS, 
with the pion beams of 10, 20, 40, 50, 80, 100, 150 and 300 GeV/c. 
The degree of description of generally accepted Bock parameteriza- 
tion of the longitudinal shower development has been investigated. 
It is shown that this parameterization does not give satisfactory de- 
scription for this combined calorimeter. Some modification of this 
parameterization, in which the e/h ratios of the compartments of 
the combined calorimeter are used, is suggested and compared with 
the experimental data. The agreement between such parameteriza- 
tion and the experimental data is demonstrated. 



1 Introduction 



One of the important questions of hadron calorimetry is the question 
of the longitudinal development of hadronic showers. This question is 
especially important for a combined calorimeter. This work is devoted to 
the study of the longitudinal hadronic shower development in the ATLAS 
barrel combined prototype calorimeter |I], |2], |3], [|, 

This work has been performed on the basis of the 1996 combined test 
beam data ||. Data were taken on the H8 beam of the CERN SPS with 
the pion beams of 10, 20, 40, 50, 80, 100, 150 and 300 GeV/c. 

2 Combined Calorimeter 

The future ATLAS experiment [||] will include in the central ( "barrel" ) 
region a calorimeter system composed of two separate units: the liquid ar- 
gon electromagnetic calorimeter (LAr) || and the tile iron-scintillating 
hadronic calorimeter (Tile) ||. For detailed understanding of perfor- 
mance of the future ATLAS combined calorimeter the combined calorime- 
ter prototype setup has been made consisting of the LAr electromag- 
netic calorimeter prototype inside the cryostat and downstream the Tile 
calorimeter prototype as shown in Fig. [l|. The dead material upstream 
of the LAr calorimeter was about 0.1 and the one between the two 
calorimeters was about 0.3 A,,-. The two calorimeters have been placed 
with their central axes at an angle to the beam of 12°. At this angle the 
two calorimeters have an active thickness of 8.6 A^. Between the active part 
of the LAr and the Tile detectors a layer of scintillator was installed, called 
the midsampler. The midsampler consists of five scintillators, 20 x 100 cm 2 
each, fastened directly to the front face of the Tile modules. The scintil- 
lator is 1 cm thick. Beam quality and geometry were monitored with a 
set of beam wire chambers BC1, BC2, BC3 and trigger hodoscopes placed 
upstream of the LAr cryostat. To detect punchthrough particles and to 
measure the effect of longitudinal leakage a "muon wall" consisting of 10 
scintillator counters (each 2 cm thick) was located behind the calorimeters 
at a distance of about 1 metre. 
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Figure 1: Test beam setup for the combined LAr and Tile calorimeters 
run. 

2.1 Electromagnetic Calorimeter 

The electromagnetic LAr calorimeter prototype consists of a stack of 
three azimuthal modules, each one spanning 9° in azimuth and extending 
over 2 m along the Z direction. The calorimeter structure is defined by 

2.2 mm thick steel-plated lead absorbers, folded to an accordion shape and 
separated by 3.8 mm gaps, filled with liquid argon. The signals are col- 
lected by Kapton electrodes located in the gaps. The calorimeter extends 
from an inner radius of 131.5 cm to an outer radius of 182.6 cm, repre- 
senting (at rj = 0) a total of 25 radiation lengths (A ), or 1.22 interaction 
lengths (A/) for protons. The calorimeter is longitudinally segmented into 
three compartments of 9 Xq, 9 Xq and 7 Xq, respectively. More details 
about this prototype can be found in Jl|, [(J . In front of the EM calorime- 
ter a presampler was mounted. The active depth of liquid argon in the 
presampler is 10 mm and the strip spacing is 3.9 mm. The cryostat has a 
cylindrical form with 2 m internal diameter, filled with liquid argon, and is 
made out of a 8 mm thick inner stainless-steel vessel, isolated by 30 cm of 
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low-density foam (Rohacell), itself protected by a 1.2 mm thick aluminum 
outer wall. 

2.2 Hadronic Calorimeter 

The hadronic Tile calorimeter is a sampling device using steel as the 
absorber and scintillating tiles as the active material pf. The innovative 
feature of the design is the orientation of the tiles which are placed in planes 
perpendicular to the Z direction f7|. For a better sampling homogeneity 
the 3 mm thick scintillators are staggered in the radial direction. The tiles 
are separated along Z by 14 mm of steel, giving a steel/scintillator volume 
ratio of 4.7. Wavelength shifting fibers (WLS) running radially collect 
light from the tiles at both of their open edges. The hadron calorimeter 
prototype consists of an azimuthal stack of five modules. Each module 
covers 27r/64 in azimuth and extends 1 m along the Z direction, such that 
the front face covers 100 x 20 cm 2 . Read-out cells are defined by grouping 
together a bundle of fibers into one photomultiplier (PMT). Each of the 
100 cells is read out by two PMTs with A0 = 2vr/64 « 0.1, while the 
segmentation along the Z axis is made by grouping fibers into read-out 
cells spanning AZ = 20 cm (At/ ~ 0.1). Each module is read out in 
four longitudinal segments (corresponding to about 1.5, 2, 2.5 and 3 A/ at 
rj = 0). More details of this prototype can be found in jj], ||. 

3 Event Selection 

We applied some cuts similar to || to eliminate the nonsingle track 
pion events, the beam halo, the events with an interaction before LAr 
calorimeter and muon events. The set of cuts applied is the following: 

• single-track pion events were selected off-line by requiring the pulse 
height of the beam scintillation counters and the energy released in 
the presampler of the electromagnetic calorimeter to be compatible 
with that of a single particle; 

• beam halo events were removed with appropriate cuts on the hori- 
zontal and vertical positions of the incoming track impact point and 
the space angle with respect to the beam axis as measured with the 
two beam chambers; 
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• a cut on the total energy rejects incoming muon. 



4 Energy Reconstruction 

To reconstruct the hadron energy in longitudinal segments the e/h 
method of the energy reconstruction, suggested in ||, has been used. In 
this method the energy of hadrons in a combined calorimeter is determined 
by the following formula: 

E = C L Ar ■ (e/ Tt)lAt ■ RLAr + CTile ' (e/ ^Tile 1 Rrile + E dm , (1) 

Here Rlat (Rnie) is the LAr (Tile) calorimeter response, clat = 1/ &LAn 
CTUe = l/criiej e LAr and e^e are the electron calibration constants for LAr 
and Tile calorimeters, e/7r ratios are equal to 

(±\ = (e/h) LAr 

^ 7T ' LAr l + ((e/h) LA r-l)L°,LAr' U 



/e \ = {e/h) Tae 

W^^e l + ((e//l)T«e-l)MT«e ' 

This method uses only the known e/h ratios and the electron calibration 
constants and does not require the previous determination of any param- 
eters by a minimization technique. The value of the (e/h) LAr ratio of 
the electromagnetic compartment has been obtained in |1(J and equal to 
(e/K)LAr = 1-77 ± 0.02. This value agrees with with estimation of > 1.7 



obtained in [ITT 



For Tile calorimeter the value of (e/K)iiie — 1.30 ±0.04 is used [12, 13 
The fraction of the shower energy going into the electromagnetic chan- 
nel for LAr compartment is 

fnO,LAr = 0.11 • ln(E beam ) . (4) 

The electromagnetic fraction in the Tile calorimeter, which samples the 
final part of shower, is equal to the one for shower with energy Exue'- 

f^\Tiie = 0.11 ■ ln(E TUe ) , (5) 

where 

E T ile = CTile ■ (e/n)Tile • Rnie ■ (6) 
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Special attention has been devoted to understanding of the energy loss 
in the dead material placed between the active part of the LAr and the 
Tile detectors. The term Ed m accounts for this. This term is taken to 
be proportional to the geometrical mean of the energy released in the last 
electromagnetic compartment (E^r^) and the first hadronic compartment 
(ETUe,i) 

Edm = Cdm ■ \J ElAt,3 ' Exile,! (J) 

similar to [|J. We used the value of Cd m = 0.31. This value has been 
obtained on the basis of the results of the Monte Carlo simulation [|TJJ]. 
These Monte Carlo (Fluka) results (open circles) are shown in Fig. ||] to- 
gether with the values (solid circles) obtained by using the expression (|7|). 
The good agreement is observed. Also the linear behaviour as a function 
of the beam energy is demonstrated. The mean energy loss is equal to 
about 3.7 ±0.4% (spread). 

The e/h method || has been tested on the basis of the 1996 test beam 
data of the ATLAS combined prototype calorimeter and demonstrated the 
correctness of the reconstruction of the mean values of energies as shown 
in Fig. |^. As can be seen the deviation from linearity for the e/h method 
is about 1%. 

We used this energy reconstruction method and obtained the energy 
depositions, E^ in each longitudinal sampling with the thickness of Axi 
in units A,,-. We transformed these depositions into the differential energy 
depositions using the formula: 

(AE/Ax)i = Ei/Axi . (8) 

Table [1| and Fig. |] show the differential energy depositions as a function 
of the longitudinal coordinate x for 10 GeV (crosses), 20 GeV (black top 
triangles), 40 GeV (open squares), 50 GeV (black squares), 80 GeV (open 
circles), 100 GeV (black circles), 150 GeV (stars), 300 GeV (black bottom 
triangles) energies. Some interesting special features are observed: the 
maximum in the region of the LAr calorimeter, then the local minimum 
in the point corresponding to the energy losses in the dead material, then 
the local maximum again. 
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5 Longitudinal Shower Development 



The next important question is understanding and description of these 
experimental data. 

There is the well known parameterization of the longitudinal hadronic 



shower development from the shower origin suggested in |15 

where Xq is the radiation length, A/ is the interaction length, a, b, d, w 
are parameters, N is the normalization factor, a = 0.6165 + 0.3193 InE, 
b = 0.2198, d = 0.9099 - 0.0237 InE, u = 0.4634. 

This parameterization is from the shower origin. But our data are from 
the calorimeter face and due to the unsufficient longitudinal segmentation 
can not be transformed to the shower origin. Therefore, we used the an- 
alytical representation of the hadronic shower longitudinal development 
from the calorimeter face []16| . This representation is a result of the inte- 
gration of the longitudinal profile from the shower origin over the shower 
position: 

d_m = UE s {x -x v ) ^ (1Q) 

(jjJb J (JjUj 



where coordinate of the shower vertex. This representation has the 

following form: 

+ (i z ^ ( . ) -^ iFi(li8+li(i _ 1) . ) } 

Here iF\(a, (3, z) is the confluent hypergeometric function and iV is the 
normalisation factor which equal to 

N= Ebeam . (12) 

A 7 r(a) (w X b~ a + (l-w) Xi d~ a ) 1 ; 

Fig. [5] demonstrates the correctness of the above formula. The calcula- 
tions by this formula are compared with the experimental data at 20 GeV 
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(crosses), 50 GeV (squares), 100 GeV (open circles), 140 GeV (triangles) 
energies for conventional iron-scintillator calorimeter |17j and at 100 GeV 
(black circles) for the Tile calorimeter as a function of the longitudinal 
coordinate x in units X^^e- The good agreement is observed. Note that 
formula (0) is given for a calorimeter characterizing by the certain X and 
A/ values. 

We have found in literature no algorithm for the hadronic shower de- 
velopment in a combined calorimeter. We suggest the following algorithm 
of combination of the LAr and Tile curves of the differential longitudi- 
nal energy deposition dE/dx (Fig. ||). The first part of the combined 
curve is the beginning of the LAr curve, the second part is the Tile 
curve. At first a hadronic shower develops in the LAr calorimeter to the 
boundary value xlav This value corresponds to certain integrated mea- 
sured energy E L a t - Then using the corresponding integrated Tile curve, 
E(x) = jQ(dE/dx)dx, (Fig. |7]) we find the point xxue in which the energy 
is equal to ETu e (xTiie) = Elat + Ed m . From this point a shower continues 
to develop in the Tile calorimeter. In principle, instead of the measured 
value of Elat one can use the calculated value of Elat = LAr (dE/dx) dx 
obtained from the integrated LAr curve (Fig. |^). 

In this way we obtained the combined curves. Fig. ^| shows a compar- 
ison between the experimental differential energy depositions at 10 GeV 
(crosses), 20 GeV (black top triangles), 40 GeV (open squares), 50 GeV 
(black squares), 80 GeV (open circles), 100 GeV (black circles), 150 GeV 
(stars), 300 GeV (black bottom triangles) as a function of the longitudinal 
coordinate x in units A„-. It can be seen that there is a significant dis- 
agreement between the experimental data and the combined curves in the 
region of the LAr calorimeter and especially at low energies. 

We attempted to understand this disagreement. We considered the 



experimental database used by Bock et al. |L5| for their parametrisation. 
It turns out that: 



The parameters given in |15| parameterization have been obtained 



by using the data from three iron-scintillator calorimeters and one 
lead-scintillator calorimeter: 

- CERN, WA1, Fe (50 mm) + Sc (6 mm), pions at 15, 50, 140 
GeV, e/h « 1.3; 

- FNAL, 379, Fe (38.4 - 51.2 mm) + Sc (6.4 mm), pions at 375 
and 400 GeV, e/h « 1.3; 
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— CERN, the combined UAl calorimeter, the electromagnetic part: 
Pb (2 - 3 mm) + Sc (1.5 - 2 mm), pions at 10, 20, 40, 60 GeV, 
e/h « 1.1 and the hadronic part: Fe (50 mm) + Sc (10 mm), 
e/ftwl.3. 

As to the used iron-scintillator calorimeters there is sufficient number 
of experimental points in our (10 - 300 GeV) beam energy range. At 
the same time the situation for the lead-scintillator calorimeter is 
quite different. There is a very limited number of points and the 
energy range of (10 - 60 GeV ) is essentially lower than used in our 
work. 

• The e/h ratio of the LAr calorimeter is ~ 1.6 times greater than 
their calorimeter. 

• This parameterization does not include such essential feature of a 
calorimeter as the e/h ratio. 

We attempted to improve the description and tried several modifica- 
tions and adjustments of some parameters of this parameterization. It 
turned out that the changes of two parameters b and w in the formula 
(pl|) in such a way b Bock = 0.22, b new = 0.34 = b Boc k • {e/h) new /(e/h) Bock , 
WBock = 0.4634, w new = 0.6 • K, where K factor is 



K=[-) /(-) , (13) 

,7T J new! \7T / Bock 



6 - = e J± (14) 

made it possible to obtain the reasonable description of the experimen- 
tal data. This is shown in Fig. || in which the experimental differential 
longitudinal energy depositions at 10 GeV (crosses), 20 GeV (black top 
triangles), 40 GeV (open squares), 50 GeV (black squares), 80 GeV (open 
circles), 100 GeV (black circles), 150 GeV (stars), 300 GeV (black bottom 
triangles) energies as a function of the longitudinal coordinate x in units 
Ajt for the combined calorimeter and the results of the description by the 
modified parameterization are compared. There is a reasonable agreement 
(probability of description is more than 5%) between the experimental 
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data and the curves taking into account uncertainties in the parametriza- 



tion function |15| . In such case the Bock parameterization is the private 
case for some fixed the e/h ratio. 

The obtained parameterization has some additional applications. For 
example, this formula may be used for an estimate of the energy deposition 
in various parts of a combined calorimeter. This demonstrates in Fig. 
|10| in which the measured and calculated relative values of the energy 
deposition in the LAr and Tile calorimeters are presented. The relative 
energy deposition in the LAr calorimeter decreases from about 50% at 10 
GeV to 30% at 300 GeV. On the contrary, the one in Tile calorimeter 
increases with the energy increasing. 



6 Conclusions 

The experimental longitudinal hadronic shower profiles in combined 
calorimeter consisting of the lead-argon electromagnetic part and iron- 
scintillator hadronic part have been obtained. The degree of description of 
the generally accepted Bock parameterization of the longitudinal shower 
development has been investigated. It is shown that this parameterization 
does not give satisfactory description for this combined calorimeter. Some 
modification of this parameterization, in which the e/h ratios of the com- 
partments of combined calorimeter are used, is suggested and compared 
with experiment. The agreement between such parameterization and the 
experimental data is demonstrated. 
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Table 1: The differential energy depositions AE/Ax as a function of the 
longitudinal coordinate x for the various beam energies. 



N 


X 


Ebearn (GeV) 


depth 


(A.) 


10 


20 


40 


50 


1 


0.294 


5.45 ±0.08 


8.58 ±0.16 


14.3 ±0.2 


16.6 ±0.4 


2 


0.681 


4.70 ±0.08 


9.10 ±0.15 


16.7 ±0.2 


20.8 ±0.3 


3 


1.026 


2.66 ±0.06 


5.55 ±0.11 


11.1 ±0.2 


13.6 ±0.2 


dm 


1.315 


1.35 ±0.07 


2.75 ±0.14 


5.28 ±0.26 


6.46 ±0.32 


4 


2.06 


1.93 ±0.03 


4.35 ±0.06 


8.99 ±0.08 


11.0 ±0.1 


5 


3.47 


0.87 ±0.02 


2.13 ±0.04 


5.29 ±0.06 


6.15±0.10 


6 


5.28 


0.18 ±0.01 


0.57 ±0.02 


1.50 ±0.03 


2.07 ±0.05 


7 


7.50 


0.025 ±0.003 


0.11 ±0.01 


0.32 ±0.01 


0.49 ±0.02 


N 


X 


E beam (GeV) 


depth 


(A.) 


80 


100 


150 


300 


1 


0.294 


22.6 ±0.6 


28.4 ±0.6 


36.3 ±0.7 


61.3 ±1.5 


2 


0.681 


30.4 ±0.4 


37.6 ±0.5 


53.5 ±0.8 


97.9 ±1.7 


3 


1.026 


20.3 ±0.3 


25.7 ±0.4 


37.2 ±0.6 


68.9 ± 1.2 


dm 


1.315 


10.1 ±0.5 


12.8 ±0.6 


19.0 ± 1.0 


34.1 ± 1.7 


4 


2.06 


18.0 ±0.1 


22.4 ±0.2 


33.9 ±0.3 


64.8 ±0.7 


5 


3.47 


11.9 ±0.1 


14.6 ±0.2 


23.3 ±0.2 


49.0 ±0.5 


6 


5.28 


3.66 ±0.06 


4.57 ±0.08 


8.18 ±0.13 


18.6 ±0.3 


7 


7.50 


0.86 ±0.03 


1.10 ±0.04 


2.04 ±0.06 


5.54 ±0.15 
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Figure 3: The demonstration of the correctness of the reconstruction of 
the mean values of energies by the e/h method (black circles): the energy 
linearity, E rec / E beam as a function of the beam energy. It is also shown the 
same for the cells weighting HI method (open circles). 
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Figure 4: The differential longitudinal energy depositions in each longi- 
tudinal sampling for 10 GeV (crosses), 20 GeV (black top triangles), 40 
GeV (open squares), 50 GeV (black squares), 80 GeV (open circles), 100 
GeV (black circles), 150 GeV (stars) and 300 GeV (black bottom trian- 
gles) energies as a function of the longitudinal coordinate x in units for 
combined calorimeter. 
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2.5 5 7.5 10 

Figure 5: The demonstration of the correctness of the formula (|ll~D- The 
comparison between the calculations by this formula and the experimental 
data at 20 GeV (crosses), 50 GeV (squares), 100 GeV (open circles), 140 
GeV (triangles) energies for conventional iron-scintillator calorimeter |17| 
and at 100 GeV (black circles) for the Tile calorimeter as a function of the 
longitudinal coordinate x in units X^^Fe 
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Figure 6: The calculated differential longitudinal energy deposition, 
dE/dx, as a function of the longitudinal coordinate x in units for the 
LAr (top) and Tile (bottom) calorimeters at 20 — 300 GeV. 
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Figure 7: The integrated longitudinal energy depositions, E(x) = 
J dE/dxdx, as a function of the longitudinal coordinate x in units of X„ 
for the LAr (top) and Tile (bottom) calorimeters at 20 - 300 GeV. 
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Figure 8: The comparison between the experimental differential energy 
depositions at 10 GeV (crosses), 20 GeV (black top triangles), 40 GeV 
(open squares), 50 GeV (black squares), 80 GeV (open circles), 100 GeV 
(black circles), 150 GeV (stars), 300 GeV (black bottom triangles) and the 
calculated curves as a function of the longitudinal coordinate x in units 

A-n-. 



19 



10 



Figure 9: The experimental differential longitudinal energy depositions 
at 10 GeV (crosses), 20 GeV (black top triangles), 40 GeV (open squares), 
50 GeV (black squares), 80 GeV (open circles), 100 GeV (black circles), 
150 GeV (stars), 300 GeV (black bottom triangles) energies as a function 
of the longitudinal coordinate x in units for the combined calorimeter 
and the results of the description by the modified parameterization. 
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Figure 10: Energy deposition (percentage) in the LAr and Tile calorime- 
ters at different beam energies. The circles (squares) are the measured 
energy depositions in the LAr (Tile) calorimeter, the diamonds (crosses) 
are the calculated energy depositions in the ones. 
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